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ABSTRACT 

We use a thin flux tube model in a rotating spherical shell of turbulent con- 
vective flows to study how active region scale flux tubes rise buoyantly from the 
bottom of the convection zone to near the solar surface. We investigate toroidal 
flux tubes at the base of the convection zone with field strengths ranging from 
15 kG to 100 kG at initial latitudes ranging from 1° to 40° with a total flux of 
10 22 Mx. We find that the dynamic evolution of the flux tube changes from con- 
vection dominated to magnetic buoyancy dominated as the initial field strength 
increases from 15 kG to 100 kG. At 100 kG, the development of fi-shaped rising 
loops is mainly controlled by the growth of the magnetic buoyancy instability. 
However, at low field strengths of 15 kG, the development of rising fi-shaped 
loops is largely controlled by convective flows, and properties of the emerging 
loops are significantly changed compared to previous results in the absence of 
convection. With convection, rise times are drastically reduced (from years to a 
few months), loops are able to emerge at low latitudes, and tilt angles of emerg- 
ing loops are consistent with Joy's Law for initial field strengths of <; 40 kG. We 
also examine other asymmetries that develop between the leading and following 
legs of the emerging loops. Taking all the results together, we find that mid- 
range field strengths of ~ 40 - 50 kG produce emerging loops that best match 
the observed properties of solar active regions. 



1. Introduction 



The Sun's cyclic large scale magnetic field with a peri od of 22 years is believed to be 
sustained by a dynamo mechanism. The Hale polarity law (jHale et al.lll919l ) of solar active 
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regions indicates the presence of a large scale subsurface toroidal magnetic field generated 
by the solar dynamo. The current prevailing picture is that the large scale toroidal magnetic 
field responsible for the formation of solar active regions is amplified and stored at the base 



of the solar convection zone (e.g. iGilmanl I200CH : ICharbonneaul 120101 ) . Thus answering the 



irp 
iuJ 



question of how active region flux tubes rise from the base of the convection zone to the 
solar surface is vitally important in the development of solar dynamo theory. 

A large body of calculations based on a so called thin flux tube approximation has 
provided important insights into the dy namic evo lution of f2-shaped rising loops in the solar 
convective envelope (see review by e.g. iFanl 120091 ) . The thin flux tube model assumes that 
the radius of the tube is significantly smaller than all other relevant scales of variation such 
as the pressure scale height and the curvature of the tube. Given the observed scale of flux 
~ 10 22 Mx in a large solar active region, and assuming that the field strength of the flux 
tubes at the base of the convection zone is at least 10 kG, which is the order of field strength 
in equipartition with convection, one finds that the above condition for the thin flux tube 
approximation is satisfied in the bulk of the solar convection zone. The thin flux tube model 
solves for the mean motion of each tube s egment und e r the influence of various i ntegrated 
forces acting on the tube segment (see e.g. ISpruitlll98ll ; lLongcope fc Klapperlll997l ). Results 
from thin flux tube models without the influence of convective flows suggest that the field 
strength of the toroidal magnetic field at the base of the solar convection zone needs to be 
in the range of about 30 kG to about 100 kG in order for the latitude of emergence and 
the tilt angles of the emerging loops to be consistent with the observed properties of solar 
active regions. If the field strength is < 20 kG , the poleward deflections of the trajectories 
of the rising flux tubes by the Coriolis force are too large such that the emer ging latitudes 
are inconsistent with the observed sunspot latitudes (e.g. ICaligari et al.lll995[ ). 



The thin flux tube calculations show that the Coriolis force acting on the rising fi-loops 
produces asymmetries between the leading and the following legs of the loop (the leading 
leg is the proceeding portion of the loop in the direction of solar rotation), which provide 
explanations for several observed asymmetries between the leading and following polarities of 
bipolar active regions. First, a slight tilt of the emerging loop is produced, with the leading 
leg being closer to the equat or than the following, consis tent with the observed Joy's law of 



active region tilt angles (e.g. iD'Silva fc Choudhurilll993[ ). Second, an asymmetric geometric 



shape of the emergi ng loop is produced, where the leading side is inclined at a smaller angle 



to the surface (e.g. ICaligari et al.l 119951 ). This may give rise to an apparent faster proper 



motion of the leading polarity as the loop emerges through the surface, which is observed 
( Ivan Driel-Gesztely fc Petrovayl ll99Cf ) . A third asymmetry is in the field strength of the 
loop, with a stronger field strength in the leading side compared to the following, which may 
explain the observed more coherent morphology of the leading polarity of an active region 
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(e.g. iFan et al.lll993f ). However, this field asymmetry is found to depen d on the assumed 
initial condition of the flux tube ( jCaligari et al.lll998l ; IFan fc Fisherlll996l ). 



The studies discussed previously do not address the dynamical effects of turbulent con- 
vection on a rising flux tube. These are shown to be adequate simplifications when the 
flux tube's magnetic field strength B > (H p / a) l l 2 B eq ~ 3 — 5B eq , where B eq ~ 10 kG is 
the field strength for which the magnetic energy density is in equipartition with the kinetic 
energy density o f convection, H p is the local pressure scale height, and a is the tube radius 



( IFan et al.ll2003l ). Solar cycle dynamo models which take into account the dynamic effects 
of the Lorentz force from the large-scale mean fields suggest that the toroidal magnetic field 

one is ~ 15 kG, and most likely cannot exceed 30 kG 



Rempel a 


2006; 


Rempel b 


2006) 



(e.g. 

flux tubes in the range of 15 kG to 30 kG rise through a turbulent solar convective envel ope- 
Thin flux tube simulations which neglect the effect of convection (e.g. see review by [Fan 
20091 ). show that tubes in this range of field strengths tend to be deflected significantly pole- 
ward by the Coriolis force during their rise, and so have difficulty reproducing the emergence 
of active regions at low latitudes, as well as Joy's law. At these field strengths, convective 
downdrafts are capable of pinning down the flux tube while helical upflows between the 
downdrafts can boost the rise of the flux tube, and therefore the resulting emerging loops 
would be significantly different from those produced in the absence of convection. 

A few studies have been performed t o investiga t e the t he buoyant rise of a fl u x tube in a 



turbu lent convective velocity field (e.g. IFan et al.l 120031 ; I Jouve &: Brunl 120091 ) . IFan et al. 



( 120031 ) carried out 3D Magnetohydrodynamic (MHD) simulations of the evolution of a 
buoyant magnetic flux tube in a stratified convective velocity field in a local Cartesian 
geometry without the effect of solar rotation. It is found that buoyant flux tubes with 
B < (Hp/aY^Beg ~ 30 kG are strongly influenced by convection, where portions of the 
tubes in convective downdrafts becomed pinned down to the bottom of the domain, while 
the rise speed of sections within upflow regions is significantly boosted. It is found that the 
evolution of the flu x tube is no longer se nsitive to the twist of the tube in this convection 
dominated regime. iJouve &: Brunl ( 120091 ) have carried out the first set of global anelastic 
MHD simulations of the buoyant rise of an initially toroidal flux ring in a rotating, fully con- 
vective spherical shell, possessing self-consistently generated mean flows such as meridional 
circulations and differential rotation, representative of the conditions of the solar convective 
envelope. Due to the limited numerical resolution and the relatively high magnetic diffusiv- 
ity in these global 3D simulations, flux tubes with a very large initial field strength (ranging 
from 150 kG to 600 kG) and a large radius, corresponding to a total flux on the order of 
a few times 10 23 Mx, significantly greater than the typical active region fluxes, are consid- 
ered. These large values are needed such that the rise times of flux tubes are less than their 
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diffusive time scales. Because most of the cases considered are essentially in the magnetic 
buoyancy dominated regime, the rising toroidal flux tube only develops rather moderate 
undulations by the influence of the convective flows, and ^-shaped tubes with undulations 
extending the depth of the convection zone are not found. 

In this study, we use the thin flux tube model, with the inclusion of a (separately com- 
puted) 3D turbulent convective velocity field in a rotating model solar convective envelope, 
to study its effects via the aerodynamic drag force, together with the forces of magnetic 
buoyancy, tension, and the Coriolis force, on the dynamic evolution of the emerging fi-loops. 
Since the thin flux tube model preserves the frozen-in magnetic field condition, it does not 
suffer the erosion of magnetic buoyancy and the tension force due to artificial diffusion of 
the magnetic field, as do the Eulerian based mult i- dimensional simulations. Although these 
thin flux tube models do not capture the internal structure and the possible fragmentation 
of the flux tube, they provide an initial step towards understanding the effects of convective 
flows in relation to other forces acting on the rising flux tubes. 

This paper is organized as follows. The equations of thin flux tube dynamics and the 
separately computed convective flow are discussed in § [2J Results of our simulations are 
presented in § [31 focusing in particular on the overall influence of convection on the flux 
tubes, the latitude of emergence and tilt angles of the flux tube, the asymmetry in field 
strength and geometry of the emerging flux loops, and the rise time for the flux tubes to 
reach the surface. We perform these thin flux tube simulations for a range of magnetic fields 
(15 kG to 100 kG) from near-equipartition to super-equipartition initial field strengths at 
varying latitudes from 1° to 40° in both the northern and southern hemispheres. A conclusion 
and discussion of the results is given in § HI 



2. Model Description 



A discussion of the thin flux tube model, and how it is derived from the MHD equations 
under the assumption that the tube radius is thin compare d to all other relevant scales of vari- 
ation can be fou nd in many previous publications (see e.g. ISpruitlll981t lLongcope fc Klapper 
19971 ; iFanl 120091 ). For this work, the equations that describe the evolution of the thin flux 
tube are: 
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where, r, v, B, p, p, T, which are functions of the time t and arc length s measured along 
the tube, denote respectively the position, velocity, magnetic field strength, gas density, 
pressure, and temperature of a Lagrangian tube segment, 1 = dr/ds is the unit vector 
tangential to the flux tube and k = d 2 r/ds 2 is the tube's curvature vector, subscript '_L' 
denotes the component perpendicular to the flux tube, $ = 10 22 Mx is the constant total 
flux of the tube, p e , p e , T e , and p, which are functions of depth only, are respectively the 
pressure, density, temperature, mean molecular weight of the surrounding external plasma, 
g is the gravitational acceleration and a function of depth, Qq is the angular velocity of the 
reference frame co-rotating with the sun, with Q set to 2.7 x rad/s in this calculation, 
Cd is the drag coefficient, set to 1 in this calculation, 7 is the ratio of specific heats, and 
v e (r,t) (discussed below) is a time dependent velocity field (relative to the rotating frame 
of reference) that impacts the dynamics of the thin flux tube through the drag force term. 
In the above equations, we do not introduce an explicit magnetic diffusion or kinematic 
viscosity term. However, we do introduce a drag force, which is the last term in Eq. [TJ 
This describes the i nteraction of th e external fluid with the flux tube in a high Reynolds 
number regime (e.g. iBatchelorl 119671 ) . The thin flux tube is untwisted, and discretized with 
800 uniformly spaced grid points along its arc length s. The numerical method used to solve 
for the flux tube evol u tion a s determined by the above set of equations has been described 
in detail in Fan et al. Jl993h . 



For the stratification of the external field free plasm a, namely p P ,, T e , /i, q, 7, and the 
super-adiabaticity, we use the reference solar model by IChristensen-Dalsgaard et al.l (119961 ) 
for the solar convection zone with an exten sion of a simple po lytropic, sub-adiabatically 



stratified thin overshoot layer, as described in lFan fc Gong] (120001 ) . Profiles of T e , p e , p e , and 



the adiabaticity 5 = V — V a d, where V = d In T e /d In p e and V a d is the value of V one 
obtains by considering local adiabatic perturbations, are shown in Figure [H The bottom 
left panel of Figure [T] shows the sub-adiabaticity of the thin overshoot layer, which extends 
from 4.8xl0 10 cm to 5.026xl0 10 cm. The bottom right panel shows the logarithm of the 
super-adiabaticity in the convection zone, which extends from 5.026 xlO 10 cm to 6.75xl0 10 
cm. 
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The main new ingredient of the current thin flux tube simulations is the inclusion into 
the drag force term of an external time-dependent convective velocity field v e (r, t) relative to 
the rotating frame of reference, computed separately from a 3-dimensional glob al convection 



simula tion using the ASH (anelastic spherical harmonic) code, as described in iMiesch et al. 



( 120061 ) (hereafter MBT06). In the anelastic approximation, convective flows and thermal 
variations are treated as a linear perturbation to a background state which is taken from 
a ID solar structure model. The computed v e (r, t) captures giant-cell convection and the 
associated mean flows (including the differential rotation and meridional circulation) in the 
rotating solar convective envelope spanning r = 0.69i? Q to r = 0.97/2© (4.8xl0 10 cm to 
6.75xl0 10 cm), resolved by a grid of 129 points in r, 256 points in 6, and 512 points in 0. 
ASH is a pseudo-spectral code with horizontal and vertical basis functions given by spherical 
harmonics and Chebyshev polynomials, eached de-aliased by keeping only the lowest 2/3 of 
modes (maximum spherical harmonic degree t max — 170 and Chebyshev degree n max = 86). 

Simulation parameters and boundary conditions are similar to Case AB3 in MBT06. 
In particular, the lower thermal boundary condition is the same as in Case AB3, with a 
latitudinal entropy gradient imposed in order to implicitly capture thermal coupling to the 
tachocline: S(9,ri) = Cp (a 2 Y" 20 + a^Y^), where S is the specific entropy per unit mass, r\ 
is the inner boundary, Cp is the specific heat at constant pressure, Ye m (6, (f>) is the spherical 
harmonic of degree I and order m, a 2 = 1.7 x 10 -6 , and = —0.43 x 10~ 6 . This helps 
promote a conical rotation profile but the strong rotational shear is attributed to the con- 
vective Reynolds stress rather than baroclinicity (see MBT06). However, the radial entropy 
gradi ent imposed at the outer bou n dary is steeper, more in line with solar structure models 



(.e.g. 



Christensen-Dalsgaard et al.l (119961 )); dS/dr = — 10 -5 erg g" 1 K _1 cm -1 in this case 
compared to — 10~ 7 erg g _1 K -1 cm -1 in Case AB3. This implies a higher Rayleigh number, 
which is partially offset by higher values of the turbulent viscosity v and thermal diffusivity 
k. In this simulation the values of v and k at the outer boundary (r = 0.97.R) are 2 x 10 13 
and 4 x 10 13 cm 2 s -1 respectively and each decreases with depth in proportion to the in- 
verse square root of the background density p -1 / 2 . The density contrast across the domain 
is 69, which corresponds to 4.2 density scale heights (somewhat larger than AB3, due to a 
slightly deeper simulation domain). This yields a mid-convection zone Rayleigh number R a 
of 5 x 10 6 and Reynolds number R e of order 50. The Rayleigh number is defined here as 
Ra = gr 2 dAS/ (ukCp) where g is the gravitational acceleration and d = r 2 — r± is the depth 
of the layer. The Reynolds number is given by V rms d/ v where V rms is the root mean square 
velocity relative to the rotating reference frame .. This value o f R a is somewhat larger than 



that used in the flux tube simulations of IJouve fc Brunl (120091 ). but R e is about a factor of 
two smaller. 



This simulation is more laminar than some others done with the ASH code (e.g. 



Miesch et al 



-7- 



20081 ; iJouve &: Brunll2009l ). and we acknowledge that this may have an impact on our results. 
However, the focus of this paper is the fundamental mechanisms by which thin flux tubes 
interact with global convection and mean flows. The convection simulation possesses all the 
relevant features necessary to investigate this interaction, including columnar, asymmetric, 
rotationally-aligned cells at low latitudes (density-stratified banana cells), a rapidly-evolving 
downflow network at higher latitudes in the upper convection zone, dominated by helical 
plumes, and a strong, solar-like differential rotation. The convective turnover timescale is 
t ~ 2xl0 6 s, or about 23 days. Our work provides a baseline for future work which can 
further assess how the details of the convective flow affect the results. Since even the highest- 
resolution simulations exhibit similar basic features, we do not expect the essential results to 
change significantly. The principle effect we expect at higher Rayleigh and Reynolds numbers 
is an increase in the random scatter due to stochastic turbulent fluctuations. Decoherence 
of large-scale motions and turbulent drag from small-scale motions could also slightly reduce 
rise times but we expect this effect to be minor since we believe that large-scale, columnar 
banana-like cells must persist even in highly turbulent parameter regimes in order to provide 
the requisite Reynolds stresses to account for the solar differential rotation. 

Figure |5] shows a snapshot of the radial velocity of the giant-cell convection at a depth of 
25 Mm below the solar surface. An associated movie showing the evolution of the convective 
flow pattern at this depth for over a period of about 315 days is also available in the electronic 
journal. The convective flow pattern shows broad upflow cells surrounded by narrow and 
intense downflow lanes. The maximum downflow speed in the convective envelope reaches 
nearly 600 m/s at a depth of about 86 Mm below the surface. Throughout most of the 
convection zone, the combined influence of the Coriolis force and the density stratification 
induces anti-cyclonic vorticity in expanding upflows and cyclonic vorticity in contracting 
downflows. This yields a mean kinetic helicity density = (u> • v) which is negative 
in the northern hemisphere and postive in the southern hemisphere (Fig. 12b), where u> 
is the vorticity of the convective flow with velocity v. There is a weak sign reversal of 
Hk in the lower convection zone where downflows expand and recirculate, inducing anti- 
cycl onic vorticity. Such a h elicity pattern is typical for rotating, compressible convection 
(e.g. 



Miesch fc Toomrel 120091 ) . 



At low latitudes there is a preferential alignment of elongated downflow lanes with the 
rotation axis, reflecting the presence of so-called "banana cells" (Fig. [2]). These features prop- 
agate in a prograde sense relative to polar regions, due in part to the d ifferential rotation 



and i n part to an intrinsic phase drift akin to traveling Rossby waves (e.g. iMiesch fc Toomre 



20091 ). Such structures dominate the convective Reynolds stress, maintaining a strong dif- 
ferential rotation comparable to that inferred from helioseismic inversions. In particular, the 
total angular velocity Q/27T (with respect to the inertia frame) decreases monotonically from 
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about 470 nHz at the equator to about 330 nHz at the poles and exhibi ts nearly conical 



conto urs at mid latitudes (see Fig. [3fr), as in the solar convection zone ([Thompson et al. 



20031 ). 



As is in iFan fc Gongl ( 120001 ). our simulations start with toroidal magnetic flux rings 
in mechanical equilibrium (neutral buoyancy), located at a radial distance to the center of 
the Sun r = r = 5.05 x 10 10 cm, slightly above the base of the solar convection zone at 
r — r czb — 5.026 x 10 10 cm. Figure H] shows a snapshot of the radial velocity at an arbitrary 
azimuthal angle 0, with the green dot representing the radius (to scale) of the largest flux 
tube, which occurs for magnetic field strengths of 15 kG. This figure shows the flux tube 
at its initial starting position compared to the base of the convection zone. Note that the 
convective velocity field is allowed to penetrate into the overshoot region. To ensure initial 
neutral buoyancy, the internal temperature of the flux tube is reduced compared to the 
external temperature. We consider a range of initial field strengths, with Bq = 15 kG, 30 
kG, 40 kG, 50 kG, 60 kG and 100 kG, and initial latitudes Ao ranging from 1° to 40° for the 
toroidal flux ring. Considering the root mean square (rms) of the convective downflows from 
the ASH simulation at the base of the convection zone which are on the order of 35 m/s, the 
equipartition magnetic field is on the order of B eq ~ 5 kG. In this case, we are investigating 
flux tubes on the order of 3 — 20 B eq . In all cases, the flux of the tube is constant at 10 22 
Mx, on the order of large scale solar active regions. For this study, we perform two sets 
of simulations sampling different time ranges of the ASH convective flow at each magnetic 
field strength for each initial latitude in both the northern and southern hemispheres. The 
toroidal ring in neutral buoyancy is perturbed with small undular motions which consist of 
a superposition of Fourier modes with azimuthal order ranging from m = through m = 8 
with random phase relations. In this work, the external time dependent 3D convective flow 
described above also impacts the flux tube through the drag force term, and we study its 
effect. 



3. Results 

3.1. Overview of Influence of Convection on Rising Flux Tubes 

Figure [5] shows snapshots of the rising flux tube in the absence of convection, developed 
from an initial toroidal flux ring at a latitude of 6°, at a time when its apex is approaching the 
top boundary, with initial magnetic field strengths of 15 kG, 40 kG, and 100 kG respectively. 
The top three images show the flux tube as if looking down on the sun from the north pole, 
whereas the bottom three images show the flux tube as if looking directly at the equator. For 
these cases in the absence of convection, emerging loops develop as a result of the non-linear 
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growth of the magnetic buoyancy instability of the initial toroidal flux tube, and the time 
for the flux tube to rise to the surface decreases with increasing magnetic field strength. The 
rise time is 6.4, 2.1, and 0.25 yrs for 15, 40 and 100 kG respectively. The rise of the loop 
becomes more radial with increasing initial field strength, with emergence latitudes of 23.3°, 
18.8°, and 7.8° for 15, 40, and 100 kG magnetic fields respectively. In Figure [5j we note the 
appearance of predominately m = 2 modes for magnetic field strengths of 40 kG and 100 
kG. The 15 kG flux tube shows a predominately m = 1 mod e, superposed with a very small 
amplitude m = 3 mode. The results are similar to those in iFan fc Gone! d200oh . where the 



toroidal flux ring develops the undular buoyancy instability (e.g. ICaligari et al.lll995l ). with 



m — 1 and m = 2 being the dominant unstable modes, and fi-shaped rising loops form. 

Figure M shows snapshots of the rising flux tubes with the same initial conditions as 
Figure but with the influence of the external convective flow included. It can be seen that 
with a low initial field strength of 15 kG, the development of the rising loops are largely 
controlled by the convective flows. Convection drastically decreases the rise time from 6.4 
yrs to only 0.2 yrs. On the other hand, at a large magnetic field strength of 100 kG, we find 
that the development of the rising loops are still mainly controlled by the development of 
the magnetic buoyancy instability, with the strongest convective downdrafts producing some 
moderate perturbations to the final emerging loops. In this case, the presence of convection 
also shortens the rise time by about one month. At 40 kG, the influence of convection is 
somewhere in between the above two cases. Again convection reduces the rise time, from 2.1 
yrs to 0.5 yrs in this case. Interestingly, with convection, the flux tube with 40 kG initial 
field takes longer to rise than either cases with 15 kG or 100 kG. The predominately m = 2 
mode for the 100 kG flux tube is still discernible. However, at lower magnetic field strengths 
the flux tube evolution is altered by convective flows such that a dominant mode can no 
longer be established. 

The differential rotation of the convective velocity field does have an influence on the 
evolution of the rising flux tubes, especially at low magnetic field strengths. Figure [7] plots 
the average velocity in the azimuthal direction of the mass elements in the legs of the 
emerging flux loop once it has reached the top of the simulation domain. We find that 
without convection, the average velocity of the mass elements of the emerging loop is always 
negative. This occurs because, as the tube rises, the conservation of angular momentum 
drives a retrograde motion of particles in the tube, producing a negative velocity in the (p 
direction. At a large magnetic field strength of 100 kG, the averaged azimuthal velocity of 
the emerging loop in the presence of convection is found to be centered around the azimuthal 
velocity for the tubes without convection. With decreasing magnetic fields, the azimuthal 
velocity in the presence of convection becomes preferentially faster than in the case without 
convection. This indicates that differential rotation is attempting to push the elements of 
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the flux tube prograde compared to the case without convection, especially at low magnetic 
field strengths where tubes are more susceptible to convective flows. 

Figure [S] summarizes how long it takes for an f2-shaped loop to emerge at the top of the 
simulation domain. Without convection, t he time for the flux tub e to emerge decreases with 



increasing field strength. This agrees with iMoreno-Insertid (119831 ). where they find the same 



trend for flux tubes given sinusoidal perturbations. Also, the range of variation in rise time 
between tubes initiating from different latitudes (as represented by the size of the error bar 
centered on the mean in Fig. [H]) decreases with increasing initial magnetic field. When the 
influence of convection is included, a different trend emerges, and the rise times are no longer 
monotonic with increasing magnetic field strength. In all cases, the rise time is reduced with 
the addition of convection, especially for the cases with weaker fields. The rise time becomes 
the longest for mid- field strength cases (i.e. 40 kG). At these mid- field strengths, the effects 
due to the magnetic buoyancy and the average convective downflows influencing the rise of 
the flux tube are of similar magnitudes, and work against each other, keeping the flux tube 
in the mid-convection zone for a longer period of time than at other magnetic field strengths. 
These effects will be discussed more in detail in the following paragraphs. For all the field 
strengths investigated, convection is found to reduce the variation of rise times for flux tubes 
with different initial latitudes, as can be seen in the size of the error bars compared to the 
case without convection. 

To understand the importance of convection, we compare the magnitude of the magnetic 
buoyancy force with that of the drag force from the external convective flows, following 



Fan et al.l ( 120031 ). For the drag force to dominate buoyancy: 

iia SnHp 

where H p is the local pressure scale height, and a is the flux tube radius a = (Q/iiB) 1 / 2 , 
with flux of the tube <J> = 10 22 Mx. The pressure scale height ranges from 5.6 xlO 9 cm at the 
base of the convection zone, to 6.9xl0 8 cm at the top of the simulation domain. Also, the 
radius of the flux tube spans from 4.6 xlO 8 cm for an initial magnetic field strength of 15 kG, 
to 1.8 x 10 s cm for 100 kG. For an order of magnitude estimate of the magnetic buoyancy 
here, we have assumed thermal equilibrium between the flux tube and the external fluid. 
Assuming 2Cd/k is of order 1, then equation (Q simplifies to 

/ \ 1/2 

V ' >V "\WJ ■ (7) 

i.e. in order for convection to dominate, the convective flow speed v e needs to be greater 
than the Alfven speed v a = B / '(47rp e ) 1,/2 multiplied by (a/H p ) 1 / 2 . 
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In Figure [HI we have plotted as a function of depth the peak downflow, peak upflow, the 
root mean square (rms) of the downflow, and the rms of the upflow of the convection velocity 
field at each constant r surface. In comparison, we have also plotted the right hand side of 
Eq. [7] evaluated at the apex (portion of flux tube with largest r value) of each of the flux 
tubes shown in Figure [6] with different initial field strengths, as it traverses the convection 
zone. It can be seen that for the 100 kG case, only the largest convective downflows are 
strong enough to impact the rise of the flux tube, whereas at 15 kG, even the rms of the 
convective flows, both upflows and downflows, can significantly affect the rise of the tube. In 
the case of 40 kG, all of the convective velocity field except for the rms upflows contribute to 
the development of the rising loops. These explain the general behavior we see in Figure [6j 
Note in Figure [9]that there is a concentration of points for the 40 kG (green) and 15 kG (red) 
cases. This is due to the fact that the flux tube is continually buffeted by convection, with 
strong downflows pushing down some apices such that new apices will take over for portions 
of the flux tube evolution. Often, the apices are also aided by boosts from convective upflows. 
Eventually, one loop becomes buoyant enough such that it will finally emerge at the top of 
the simulation domain. There is a jump in the 15 kG curve that occurs near 6.3 xlO 10 cm 
because one loop is boosted by a convective upflow and becomes significantly more buoyant 
than the previous apex. It is clear from Figure M, that downflows in the convective velocity 
field dominate in amplitude, however their spatial extent is small compared to the upflows as 
indicated by the narrow downflow lanes shown in Figure [21 Strong downflows can pin the flux 
tubes to the base of the convection zone at the beginning of the simulation faster than the 
Fourier mode perturbations can in the case without convection. Especially at low magnetic 
field strengths where the tube is highly susceptible to convection deformation, many portions 
of the tube will become anchored as compared to the case without convection. Also, rising 
loops can be significantly boosted by broad upflows such that they emerge at the top of the 
simulation domain much faster than in the case without convection, where only buoyancy 
aids in driving the tube to the surface. Convection can also enhance buoyancy instabilities 
by introducing finite-amplitude perturbations and subsequent gravity induced draining of 
fluid from the flux tube apex. Figure [10] shows snapshots (one individual time instance) of 
the radial distance r of the flux tube from Sun center as a function of the azimuthal angle 
(f) (black lines), as well as the external radial velocity experienced by the tube at the height 
r of each tube segment (red lines). These snapshots are for the same flux tubes as used in 
Figures [6] and [9j with the left column showing the snapshot at a time when the apex of the 
tube has reached O.82i? , and the right column showing the snapshot for the last time step 
when the flux tube apex has reached the top of the simulation domain at 0.97i? Q . From 
this figure, it is evident that at large magnetic field strengths, only the strongest downflows 
can perturb the tube. However, at small magnetic field strengths, all flows are capable of 
deforming the flux tube. 
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3.2. Latitude of emergence and tilt angles 

In the absence of convection, we find that for weaker initial field strengths, there is a low 
latitude zone devoi d of flux emergence due to the poleward deflectio n of the rising tube by the 



Coriolis force (e.g. IChoudhuri &: Gilmanlll987t ICaligari et al.lll995f ). This is shown in Figure 



[TT1 where we have plotted the deflection of the flux tube in latitude (emergence latitude 
minus initial latitude) as a function of the flux tube initial latitude. The red diamonds (blue 
plus signs) show the deflection of the flux tubes without (with) the influence of convection. 
This low latitude zone void of flux emergence in the case without convection shrinks as the 
initial magnetic field strength is increased. This is because at larger initial magnetic field 
strengths, the buoyancy force overpowers the Coriolis force such that tubes rise more radially. 
However, when the flux tubes are subjected to the convective flow, we find that they are 
able to emerge at lower latitudes near the equator, even for low field strength cases. In some 
cases, the convective velocity field is able to push the flux tube such that its apex will emerge 
closer to the equator than where the tube initially started. This is indicated in Figure [TT] 
by negative plus sign points. We note that with convection, the poleward deflection of the 
flux tubes is reduced for the majority of cases. Flux tubes subjected to turbulent convection 
spend less time in the convection zone than the flux tubes without convection (see Fig. [8]). 
As such, the Coriolis force will have less time to deflect these flux tubes to higher latitudes. 
Also, increased anchoring of weaker flux tubes by convection reduces the fraction of the tube 
that moves away from the rotation axes, and therefore reduces the effect of the Coriolis force 
on the tube. 

In most cases without convection, the variation of the emergence latitude with initial 
latitude is smooth, except for the 15 kG flux tube, where there is an abrupt 14° jump in 
emergence latitude between the tubes originating at 5° and 10° (see diamond points in the 
bottom right panel of Fig. [11] and Fig. fT2l) . At this magnetic field strength in the absence 
of convection, flux tubes originating at 5° and 1° have troughs which become anchored to 
the overshoot region below the base of the convection zone. However for the same tubes 
originating from 8° and above, the rising tubes do not become anchored. As such, the tube 
floats upward as a whole, and experiences a relatively greater deflection towards the pole 
and hence emerges at a higher latitude compared to the anchored cases. In these cases, the 
m = 1 to m = 3 modes grow slower than the m = mode. This corresponds to a poleward 
slip of the tube as a whole. Thereafter, modes on the order of m = 1 to m = 3 begin to 
develop larger amplitudes, allowing loops to develop, though the troughs do not anchor at 
the base of the convection zone. With convection, at this same magnetic field strength, all 
of the flux tubes have troughs which become anchored because convective downflows pin 
them to the bottom of the convection zone from the beginning of the simulation. We note 
that the flux tubes tend to anchor (when they anchor at all) around 5.0xl0 10 cm in both 
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the cases with and without convection. At this depth, the sub-adiabaticity is 5 =-1.3 xlO -4 . 
The flux tubes all originate at 5.05 xlO 10 cm, 2.8 xlO 8 cm above the base of the convection 
zone, and have a super-adiabaticity here of 5 =3.2 xlO -8 , or very nearly zero. Reducing 
the sub-adiabaticity would reduce the chances for anchoring of the portions of the flux tube 
that are perturbed into the overshoot region. Increasing the sub-adiabaticity would help 
anchor these flux tubes if they were perturbed into the overshoot region and did not remain 
anchored at a reduced sub-adiabaticity. However, we find that the tubes which do not anchor 
in this study are never perturbed into the overshoot region. Instabilities of higher order with 
large enough amplitude to perturb the tube into the overshoot region do not develop fast 
enough to overcome the m = mode. 

Figure [T2] shows the tilt angle of the emerging flux loops as a function of the emergence 
latitude, for flux tubes with a range of initial latitudes (from 1° to 40° as indicated by the 
color of the points). The different panels in the Figure show the results for initial magnetic 
field strengths of 15 kG, 30 kG, 40 kG, 50 kG, 60 kG, and 100 kG respectively. In the plots, 
results with (without) the influence of convection are shown as plus signs (diamonds). 

It is well known that bipolar active regions on the sun tend to form with the leading 
polarity (in the direction of solar rotation) being slightly closer to the equator than the 
following polarity, and thus showing a tilt angle for the line connecting the following to the 
leading polaritie s with respect to the eas t west direction, as described by Joy's Law (e.g. 



Hale et al.l Il919l ; iDasi-Espuig et al.ll2010l ). On average, the tilt angle of the active region 
increases with its latitude. In Figure [T2] the tilt angle is computed as the angle between 
the tangent vector at the apex of the emerging loop and the local east-west direction. We 
define a positive sign of tilt as a clockwise (anti-clockwise) rotation of the tangent vector 
away from the east-west direction in the northern (southern) hemisphere, consistent with 
the direction of the observed mean tilt of active regions. If on the other hand the magnitude 
of the tilt angle exceeds 90°, then it means that there is an anti-Hale polarity arrangement 
for the active region. 

It is important to note that we do not consider the effect of magnetic field line twist 



on the evolution of flux tubes. iFanl (120081 ) considers the effect of twist on the tilt angle 
of emerging flux tubes for a 3D MHD simulation in a local Cartesian geometry without 
an external velocity field. They find that in order for emerging flux tubes of highly super- 
equipartition field strength on the order of 10 5 G to have tilts consistent with Joy's Law, 
the initial twist rate of the flux tube needs to be smaller than about a half of that required 
for the tube to rise cohesively. However, in the case of flux tubes in a convecting box with 
equipartition magnetic field strengths on the order of 10 4 G, it is found that convective 
downflows dominate over magnetic buoyancy effects, and the evolution of the tube is no 
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longer sensitive to the initial twist of the tube ( IFan et all 120031 ). In this paper, we are 
primarily concerned with how convection alters the overall evolution of the rising flux tubes, 
and the thin flux tube model used here does not address the effect of the field line twist. 



Dasi-Espuig et all (120101 ) analyzed Mount Wilson and Kodaikanal sunspot group tilt 
angle and latitude data spanning solar cycles 15-21. In order to re-derive Joy's Law, they 
performed a fit for the equation: a = mX to the data, where a, m, and A represent the tilt 
angle, the slope, and the latitude respectively. The fit is forced to go through zero because 
no tilt is expected for equatorial sunspot groups. They obtained a slope of 0.26 ± 0.05 and 
0.28 ± 0.06 for Mount Wilson and Kodaikanal data, respectively. 

In Figure [121 we show the same kind of linear least squares fit for both the simulated tilt 
angles with and without the influence of the convective flow. The slope of the best fit lines 
with uncertainties for 100 kG, 60 kG, 50 kG, 40 kG, 30 kG, and 15 kG without convection are 
respectively: 0.25±0.01, 0.30±0.02, 0.22±0.02, 0.12±0.01, -0.0 7±0.05, and -0.13±0.i3 . 



Thus the slopes of the best fit lines fall within the range found by iDasi-Espuig et all (120101 ) 
only for flux tubes with initial field strengths of 100 kG, 60 kG, and 50 kG. For weaker 
initial field strengths (below about 40 kG) the emerging loops begin to show negative tilt 
angles, opposite to the sign of the active region mean tilts. This occurs because flow along 
the flux tube at the loop apex cha nges from diverging into converging as i t enters the upper 



layers of the convection zone (e.g. ICaligari et al.lll995l ; IFan fc Fisherlll996l ). and the Coriolis 



force acting on the converging flow drives a tilt of the wrong sign. With the influence of 
the convective flow included, the slopes of the best linear fits to the tilt angles all increase 
due to the systematic effect of the kinetic helicity (see Fig. [3]^) associated with the rotating 
convective flows, which tends to drive rotation at the apex of the loops in the direction 
consistent with the sign of the observed active region mean tilts. The slope of the best fit 
lines with convection for initial magnetic field strengths of 100 kG, 60 kG, 50 kG, 40 kG, 30 
kG, and 15 kG are respectively: 0.30 ±0.08, 0.34 ±0.05, 0.34 ±0.12, 0.44 ±0.14, 0.12 ±0.18, 
and 0.15 ± 0.21. As a result, the tilt angles of emerging loops with initial field strengths of 
100 kG, 60 kG, 50 kG, and 40 kG are all consistent with Joy's Law. These field strengths 
all show a positive Joy's Law trend, given the uncertainties in the fitted slopes. However, at 
lower field strengths, the tilt angles show large random scatters produced by convection, and 
do not show a significant systematic dependence on latitude as described by Joy's law. The 
uncertainty of the slope also tends to increase with decreasing magnetic field, reflecting the 
random scatter produced by convection. While the slope of the best fit lines for 30 kG and 
15 kG are still positive, the uncertainties are too large to report a definitive positive Joy's 
Law trend for these field strengths. There is a larger than expected uncertainty on the slope 
for 100 kG. This occurs because there is an outlier, which can not be seen in Figure [12] as it 
falls outside the region graphed. The outlier corresponds to a flux tube which originates in 
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the southern hemisphere, but emerges in the northern hemisphere such that the emergence 
of the tube would result in an anti-Hale polarity arrangement for the active region. 



3.3. Asymmetry in Field Strength 



Another well-known observed asymmetry is in the morphology of the leading and the 
following polarities of an active region, where the leading polarity flux tends to be concen- 
trated into a well formed sunsp ot, whereas the following p olarity flux tends to appear more 



fr agmented and dis persed (e.g. iBray fc Loughheadl Il979l ). The thin flux tube simulations 



of iFan et al.l (119931 ) showed that the preceding leg of an emerging magnetic flux loop has 
a stronger magnetic field than the following leg as a result of the differential stretching of 
the rising loop due to the Coriolis force. This gives an explanation for the observed more 
coherent and less fragmented morph ology for the leading polarity flux in an active region. 



However subsequent simulations (e.g. ICaligari et al.lll995l : IFan fc Fisherlll996l ) using the me- 



chanical equilibrium i nitial state (whic h is more physical) as opposed to the temperature 
equilibrium as used by Fan et al.l ( 119931 ) , found that the leading leg of the emerging loop has 
a stronger magnetic field than the following only for flux tubes with an initial field strength 
that is below about 60 kG. For flux tubes with higher initial field strengths, the field strength 
asymmetry reverses at the top of the loop. 

Here we investigate this magnetic field asymmetry by calculating dB/ds at the apex of 
the emerging flux loop as shown in Figure [TBI If dB/ds is greater (less) than zero, then the 
leading (following) leg has a stronger magnetic field. As can be seen in Figure HU at 50 kG 
and below, the majority of the emerging loops show stronger field in the leading leg than 
the following leg in the presence of convection. However, at 100 kG, with convection, only 
some of the emerging flux loops at higher latitudes have stronger field in the preceding leg. 
At 60 kG there are about equal number of positive and negative dB/ds cases. Thus with 
convection, loops with initial field < 50 kG tend to emerge with the appropriate magnetic 
field asymmetry (i.e. with a stronger field in the preceding leg and therefore can be expected 
to result in an emerging active region with a more coherent leading polarity). 



3.4. Asymmetry of Inclination 

When looking down on the emerging loops from the pole, it is apparent that there is an 
asymmetry in the inclination of the leading and following legs with respect to the vertical 
direction, as is evident in the cases without convection shown in Figure The following leg 
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tends to have a steeper slope than the leading leg. This asymmetry is caused by the Coriolis 
force and the conservation of angular momentum as the tub e rises through the convection 
zone (IMoreno-Insertis et al.lll994t ICaligari et al.l Il995l . Il998l ). and provides an explanation 
for the apparent asymmetric east-west proper motions of the two polarities of an emerging 
active region. As the tube rises above a constant r surface, inclination of the legs of the loop 
cause an apparent more rapid motion of the leading polarity spot away from the emerging 
region as compared to the motion of the following polarity spots. 

Here we quantify the steepness of each leg of the emerging loop by calculating what we 
call the inclination angle. To do this, first we find the portion of the emerging loop which 
is concave downward. Then we find the best fit line for each leg from the flux tube apex 
in the concave down portion. The angle between the best fit line for each leg and — f is 
the inclination angle. The smaller this angle, the steeper the slope of the leg. In Figure 
HH we have plotted the average difference of the inclination angles (leading leg minus the 
following leg). Only the average inclination difference for each magnetic field strength is 
plotted because we find that the inclination difference does not vary systematically with 
latitude. The bars are the standard deviation, showing the spread in inclination angle 
difference for each magnetic field strength. This spread is due primarily to convective effects. 
Po sitive inclination difference m eans the following l e g is s t eeper , consistent with the results 
of IMoreno-Insertis et all rtl994T ) and ICaligari et all jl995l Il998h . We find that at all field 
strengths, the majority of the emerging loops develop a steeper slope for the following leg. 
The inclination angle differences are overall reduced with the inclusion of convection. 

It is worthwhile to note that convection aids in bringing the footpoints of the rising flux 
loops closer together, most notably in the 15 kG and 40 kG cases as can be seen in Figure 
[5] and Figure |6j This can also be observed in Figure [TUl where strong convective downflows 
pin portions of the tube down, creating many buoyant loops with very close footpoints at 15 
kG. 



4. Discussion 

In subjecting the thin flux tube model to a turbulent solar convective velocity field 
computed separately from a 3D global convection simulation in a rotating spherical shell 
representing the solar convective envelope, we study how convection can influence the de- 
velopment and emergence of solar active region flux tubes. Although idealized, the thin 
flux tube approximation allows us to investigate flux emergence at moderate to strong field 
strengths (15 - 100 kG) in localized flux concentrations under perfect flux frozen-in condi- 
tions. We find that as the field strength of the initial toroidal flux ring placed near the 
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bottom of the solar convective envelop is increased from 15 kG to 100 kG, the develop- 
ment and the evolution of the emerging flux loops change from being convection dominated 
to magnetic buoyancy dominated. At 15 kG initial field strength, the development of the 
emerging flux loops is largely controlled by the convective flows. Both the average up-flows 
and down-flows can significantly affect the rise of the flux tube. On the other hand, at 100 
kG initial field strength, the development of the emerging loops are largely controlled by the 
growth of the magnetic buoyancy instability and only the strongest convective downdrafts 
can significantly impact the rising tubes. 

With the inclusion of the convective flows, the previous issue of 15 kG - 30 kG flux tubes 
being significantly deflected poleward by the Coriolis force during their rise is resolved. Loops 
can now emerge at low latitudes all the way to the equator. The main reasons for this are the 
anchoring of the emerging flux loops by the convective downdrafts and also the faster rise of 
the loops propelled by the upflows. Convection is also found to produce random scatters in 
the tilts of the emerging loops, especially for the weaker field strength cases where the scatters 
are greater. In addition, because the convective flow in the rotating spherical envelope shows 
a mean kinetic helicity that is negative (positive) in the northern (southern) hemisphere, it 
on average tends to drive clockwise (anti-clockwise) tilts for the rising loops in the northern 
(southern) hemisphere, consistent with the sign of the active region mean tilts. As a result, 
the inclusion of convection tends to increase the slope of the linear Joy's Law fit of the tilt 
angles of the emerging loops as a function of emerging latitudes. Flux tubes with initial field 
> 40 kG and < 100 kG, are found to produce a mean tilt angle dependence on latitude 
that is consistent with the observed Joy's Law. However the effect of the kinetic helicity in 
the convective flow is not sufficient, on average, to correct the tendency for loops to develop 
tilts of the wrong sign (opposite to that of the observed active region mean tilt) for tubes 
with initial field strength < 40 kG. 



Similar to ICaligari et all ( 119951 ) and Fan fc Fisherl ( 119961 ). we find that at or below 



50 kG, the leading leg of the emerging loop tends to have a larger magnetic field than 
the following, which may provide an explanation for the observed better cohesion of the 
leading polarity of an emerging active region compared to the following polarity. This trend 
of asymmetry in field strength reverses for tubes with initial field > 60kG. Also we find 
that the inclusion of convection does not change the general trend of the asymm etry in the 



geometry of the emerging loop found in previous studies (e.g. ICaligari et al.lll9951 ). where the 
leading leg tends to be more horizontally inclined than the following, which would provide 
an explanation for the apparent asymmetric east-west proper motions of the two polarities 
of an emerging active region. This trend of geometric asymmetry is not dependent on the 
magnetic field strength. 
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These results combined suggest that in order for the emerging loops to be consistent with 
the observed properties of solar active regions, the best range of field strengths for toroidal 
magnetic field at the bottom of the convection zone is ~ 40 - 50 kG. Above this range, the 
magnetic field is strong enough to resist convection such that the flux tube properties are not 
changed significantly from the case without convection, but the magnetic field asymmetry of 
the emerging loops tend to be of the wrong sense. Below this range, convection dominates 
the development and evolution of the emerging loops. The loops are able to emerge at low 
latitudes, but their tilt angles show very large scatter. Increasing the number of simulations 
run to collect more data may result in a reduction of the uncertainty on the slope of the 
Joy's Law best fit lines, such that 15 kG and 30 kG flux tubes may also show a positive 
Joy's Law trend. 

The main limitation of the thin flux tube approach is that the model assumes the 
tube remains cohesive, and does not address the dynamic effect of the internal twist of 
the flux tube. Also, as the ASH convective flow is computed separately, the calculations 
do not address how the rising flux tube might affect the surr ounding convective flow. Full 



3D spherical-shell MHD simulations (e.g. iJouve fc Brunl 120091 ) are necessary to address the 



above questions, but it is still beyond the current computational capabilities to adequately 
resolve an active region scale flux tube in a global convection simulation. It is important to 
note that these results are for flux tubes which have reached 21 Mm below the solar surface. 
Even though these flux tubes still must traverse a small portion of the convection zone and 
encounter the solar surface shear layer, we feel that the behavior of the flux tubes at this 
depth may give a good representation of the large scale pattern of flux emergence at the 
solar surface. This paper serves as an initial study to understand how active region scale 
flux tubes of various strengths are affected by convective flows as they rise through the solar 
convective envelope. 
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Fig. 1. — Profiles of T e , p e , p e and adiabaticity factor 5 for the entire simulation domain as 
a function of solar radius r for the reference solar model. The bottom left panel shows the 
sub-adiabaticity in the overshoot region, whereas the bottom right panel shows the logarithm 
of the super-adiabaticity in the convection zone. 
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Fig. 2. — Snapshot of convective radial velocity at a depth of 25 Mm below the solar surface. 
A movie showing the evolution of the radial velocity at this depth for over a period of about 
315 days is available in the electronic version of the paper 
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Fig. 3. — (a) Mean kinetic helicity density and (b) angular velocity (with respect to the 
inertial frame) in the convection simulation, averaged over longitude and time (time interval 
1366 days). Color tables saturate at the values indicated, with extrema ranging from (a) 
-0.122-0.133 cm s^ 1 and (b) 326-468 nHz. 



Fig. 4. — Snapshot of a slice of the radial velocity field (r = 0.69i? Q to O.97i? , 9 from north 
solar pole to equator) at an arbitrary azimuthal angle 0, with a cross-section of a 15 kG flux 
tube (green, drawn to scale), at its initial starting radius 6° above the equator. The dotted 
line represents the base of the convection zone at 5.026 xlO 10 cm. 
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B= 15 kG B= 40 kG B= 100 kG 

t (days) - 2321.4 t (days) - 778.2 t (days) = 91.4 




B= 15 kG B= 40 kG B= 100 kG 

t (days) = 2321.4 t (days) - 778.2 t (days) = 91.4 



Fig. 5. — Snapshots of the rising flux tube, developed from an initial toroidal flux ring at a 
latitude of 6°, at a time when its apex is approaching the top boundary, with initial magnetic 
field strengths of 15 kG, 40 kG, and 100 kG respectively, in the absence of convection. The 
top images show a polar view, whereas the bottom images show an equatorial view. In all 
cases, the image has been rotated such that the flux tube apex is on the right, and at the 3 
o'clock position if looking down from the north solar pole. The orange sphere has a radius of 
4.9 xlO 10 cm, which is O.lxlO 10 cm above the base of the rotating solar convective envelope, 
and 0.15xl0 10 cm below the initial position of the flux tubes. 



B= 15 kG B= 40 kG B= 100 kG 

t (days) = 73.1 t (days) - 195.4 t (days) - 62.4 




B= 15 kG B= 40 kG B= 100 kG 

t (days) = 73.1 t (days) - 195.4 t (days) 62.4 



Fig. 6. — Same as Figure El except the tube is subjected to the external convective flow. 
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Fig. 7. — Average velocity in the azimuthal direction <p of the mass elements in the legs of 
the emerging flux loop once it has reached the top of the simulation domain. These values 
are plotted for initial magnetic field strengths of lOOkG, 60 kG, 50 kG, 40 kG, 30 kG, and 
15 kG, respectively, in the presence of convection (plus signs), and the absence of convection 
(diamonds). 
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Fig. 8. — Average time for the flux tube to rise through the convection zone for magnetic 
fields of 100 kG, 60 kG, 50 kG, 40 kG, 30 kG, and 15 kG in the absence of convection (red 
diamond points), and in the presence of convection (blue star points). Bars represent the 
variation in the rise time for different starting latitudes of the flux tube. 
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Fig. 9. — The peak down-flow, peak up-flow, the root mean square (rms) of the down-flow, 
and the rms of the up-flow of the convection velocity field in each constant r surface as a 
function of r. Also plotted are the right hand side of Eq. [7] at the apex of each of the the flux 
tubes shown in Figure [6] with different initial field strengths, as it traverses the convection 
zone. 
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Fig. 10. — Flux tube radial distance from center of the Sun r (black line), plotted with the 
external radial velocity experienced by the flux tube at the height r of the flux tube segment 
(red line), both as functions of the azimuthal angle <fi. Figures on the left are for an instance 
when the apex height is r = 0.82i? Q , and on the right are for the last time step when the 
apex has reached the top of the simulation domain at r = O.97_R . These snapshots show 
all the azimuthal angle values for the flux tube, with all tubes originating at 6° latitude. 
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Fig. 11. — Latitudinal deflection (emergence latitude minus initial latitude) of flux tube 
apex as a function of initial latitude for initial tube field strengths of 100 kG, 60 kG, 50 
kG, 40 kG, 30 kG, and 15 kG in the absence of convection (red diamond points), and in the 
presence of convection (blue plus sign points). 
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Fig. 12. — Tilt angle as a function of emergence latitude for initial magnetic field strengths of 
100 kG, 60 kG, 50 kG, 40 kG, 30 kG, and 15 kG respectively, and for cases with (plus signs) 
and without (diamond points) the influence of convection. The gray line is the least squares 
fit for flux tubes in the absence of convection. The best fit line slopes with uncertainties for 
field strengths of 100 kG, 60 kG, 50 kG, 40 kG, 30 kG, and 15 kG are 0.25 ±0.01, 0.30 ±0.02, 
0.22 ±0.02, 0.12 ±0.01, -0.07 ±0.05, and -0.13 ±0.13. The black line is the fit to the cases 
subjected to the convective flow, with best fit line slopes for fields strengths of 100 kG, 60 
kG, 50 kG, 40 kG, 30 kG, and 15 kG of 0.30 ± 0.08, 0.34 ± 0.05, 0.34 ± 0.12, 0.44 ± 0.14, 
0.12 ± 0.18, and 0.15 ± 0.21. A color bar indicates the original starting latitude of the flux 
tube. Dashed lines indicate 0°. 
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Fig. 13. — dB/ds at the apex of the emerging flux loop as a function of emergence latitude 
for tubes with an initial field strength of 100 kG, 60 kG, 50 kG, 40 kG, 30 kG, and 15 
kG respective, for cases with (plus signs) and without (diamond points) the influence of 
convection. A color bar indicates the original starting latitude of the flux tube. Note that 
the y-axes are not the same for every plot. 
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Fig. 14. — Average difference of the inclination angles (see description in the text) between 
the leading and the following sides of the emerging loop for initial magnetic field strengths 
of 100 kG, 60 kG, 50 kG, 40 kG, 30 kG, and 15 kG respectively, for cases with (blue star 
points) and without (red diamond points) the influence of convection. Bars represent the 
variation in inclination angle difference for different starting latitudes of the flux tube. 



